Draft version July 14, 2008 

Preprint typeset using L^T^X style cmulatcapj v. 03/07/07 



00 
O 
O 
(N 



SPITZER AND HHT OBSERVATIONS OF BOK GLOBULE B335: ISOLATED STAR FORMATION 

EFFICIENCY AND CLOUD STRUCTURE l 

Amelia M. Stutz 2 , Mark Rubin 3 , Michael W. Werner 3 , George H. Rieke 2 , John H. Bieging 2 , Jocelyn Keene 4 
Muu Kang 2,5,6 , Yancy L. Shirley 2 , K. Y. L. Su 2 , Thangasamy Velusamy 3 , David J. Wilner 7 

Draft version July 14, 2008 

ABSTRACT 

We present infrared and millimeter observations of Barnard 335, the prototypical isolated Bok 
globule with an embedded protostar. Using Spitzer data we measure the source luminosity accurately; 
we also constrain the density profile of the innermost globule material near the protostar using the 
observation of an 8.0 (im shadow. HHT observations of 12 CO 2-1 confirm the detection of a flattened 
molecular core with diameter ~ 10000 AU and the same orientation as the circumstellar disk (~ 100 to 
200 AU in diameter). This structure is probably the same as that generating the 8.0 /xm shadow and 
is expected from theoretical simulations of collapsing embedded protostars. We estimate the mass of 
the protostar to be only ~5% of the mass of the parent globule. 
Subject headings: ISM: globules - ISM: individual (Barnard 335) - infrared: ISM - stars: formation 
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iBok fc Reillvl (J1947D and [Bok| (|1948ft were the first to 
suggest that "Bok globules" could be the sites of iso- 
lated star form ation. With this in m i nd, B 335 (CB199 
in the catalog of lClemens fe Barvainis! (11988D1 was exten- 
sively studied i n molecular lines (jMartin fc Barrct tjfl978l: 
Ipickman1ll978f) an d in the f ar-infrared bv | Keene et all 
(|1980T ) and iKeend rfl98l . iFrerking fe Langeri (11982ft 
discovered a "pedestal" source, i.e., a molecular out- 
flow, indicating the possible presence of a protostar 
within the cloud. The presence of a protostar was 
quickly co nfirmed with addit ional far-infrared obser- 
vations by iKeene et al.l (|1983l). Assuming a dis tance 
of 250 pc, as suggested by iTomita et al.1 (Il979f) and 



as ado pted by most subsequent observers, IKeene et al.l 
(1983) determined the luminosity of the source within 
B335 to be about 3 L . B335 has subsequently 
been observed at far-infrar e d and submillim eter wave- 
length s by, e.g., iGee et ail (119851). iMozurkewich et al.1 
(fl98ffl.Oi andler e t all (J1990D . iHuard et al.l (119990 . and 
iShirley et al.l (j2000f k it is the first example of a Class 
young stellar object to be so thoroughly studied. It is 
therefore the template source for many of our concepts 
about the early stages of the formation of low-mass stars. 
Modeling has established that the globule material 
is still undergoing str ong noncircular motions that are 
most likely infall (e.g IZhou et al.l[i~990l : IChoi et al.lfl995l : 
IVelusamv et al.lll995[ ). Although this scenario of infall, 
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and more specifically of inside-out collap se, is plausi- 
ble, important discrepancies remain (e.g.. JWilner et al.l 
[20001: lHarvev et al1l200ll : lEvans et al.ll2005l : lChoill2007ir 

Models and obser vations imply tha t the core is a star 
of mass ~ 0.5 M Q (Zhou et al. 1990) with a circumstel- 
lar dis k ~ 100 to 200 AU in diameter (|Harvev et al.l 
2003al). This source ha s a strong bipolar ou tflow (e.g., 



Cabrit fc BertoutJ 119921 : iHirano et all [1988), oriented 



nearly in the plane of the sky, that maintains a cavity in 
the interstellar medium (IS M) and is associated with the 
observed shocked material (jGalfalk fc O lofsson 2007). 

We describe observations with all three instruments 
on Spitzer and with the Heinrich Hertz Submillimeter 
Telescope that allow us to 1) measure the density dis- 
tribution in the center of the globule; 2) constrain the 
geometry of the outflow and thus the protostellar disk; 
and 3) measure the source luminosity accurately. These 
measurements (with a revised distance of 150 pc) reveal a 
flattened molecular cloud core, consistent with numerical 
simulations of cloud collapse; provide tight constraints 
on the opening angle of the outflow; and demonstrate in- 
ternal consistency between the globule luminosity (and 
other characteristics) and those expected for a 0.3 to 
0.4 M Q young star. 

2. OBSERVATIONS AND PROCESSING 

2.1. Spitzer data 

B335 was observed by the Spitzer Space Telescope on 
2004 April 20 and 21 with the Infrared Array Cam- 
era fIRAC: iFazio et al.l l200l AORKEYs 4926208 and 
4926464; 2004 October 15 with the Multi band Imag- 
ing Photometer (MIPS; iRieke et all l200l AORKEY 
12022106; and 2004 October 14 w ith the Infrared Spec- 
trograph (IRS; iHouck et all l200l AORKEY 3567360. 
Observations were made at two closely separated epochs 
to permit the elimination of asteroids from the final mo- 
saic. The field of view with coverage in all four IRAC 
bands is 5.5' x 13.6', with a larger field covered in at least 
two bands. For each IRAC band, 3 overlapping positions 
arranged along a column were imaged. At each posi- 
tion five small Gaussian dithers of 10.4 s duration were 
taken. In the region where we had coverage in all four 
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Fig. 1. — Gallery of images of B335 (CB199) at the indicated wavelengths. Images are shown with a logarithmic intensity scale and are 
1'6 on a side. North is up and east is to the left. 



bands, this resulted in an average exposure time of about 
120 sec. For each dither, the 10.4 second long-integration 
frame was accompanied by a short-integration frame of 
0.4s in "high dynamic range (HDR)" mode. This allowed 
for correction of saturated pixels in the longer exposures. 

MIPS instrument observations were carried out using 
scan map mode in all three MIPS bands at 24 /zm, 70 ^m, 
and 160 /xm as part of Spitzer program ID 53 (P.I. G. 
Rieke). The data reductio n was carried out u sing the 
Data Analysis Tool (DAT; iGordon et aT1l2005h o ut fol- 
lowing the steps outlined in a paper bv lStutz et al.l (2007) 
describing a similar study of the Bok Globule CB190; 
please refer to that publication for details. For B335, 
the observations covered a held of ~ 15' x 55', with an 
exposure time of ~ 192 s at 24 fim, ^80 s at 70 /<m, and 
~ 20 s at 160 /xm, although the integration time varies 
significantly in each map due to non-uniform coverage. 

The IRAC data were processed by the Spitzer Science 
Center (SSC) using the standard data pipeline (version 
11.4) to produce basic calibrated data (BCD) images. 
To process the BCDs and produce mosaics, we used a 
set of IDL procedures and csh scripts (for data handling 
and module control) that were provided by Sean Carey of 
the Spitzer instrument support team. These routines use 
existing MOPEX software developed by the Spitzer Sci- 
ence Center (Mak ovoz fc Marleaull2005r ). We produced 
final mosaics with a pixel scale of 076 (see Figure [1]). 
Source extraction and ph otometry were performed using 
the SExtractor software (|Bertin fc Arnoutsl [T996). An 
aperture with a diameter of 8 pixels was used in the 
extraction. Aperture correction factors, derived by the 
Spitzer Science Center and provided in the IRAC Data 
Handbook 2.0, were used to estimate the true flux of ob- 



jects based on the measured flux in each aperture. These 
fluxes were then converted into magnitudes using the zero 
magnitude flux values taken from the IRAC Data Hand- 
book. The overall flux calibration, expressed in physical 
units of MJy sr _1 , is estimated by the SSC to be accu- 
rate to within 10%. All the Spitzer images are displayed 
in Figure [TJ even at 8.0 /xm the embeded point source is 
not visible and the images are dominated by reflection 
nebulosity. The Spitzer photometry is plotted in Fig- 
ure [2] along with longer wavelength data taken from the 
literature. 

IRS observations of B335 were taken in mapping mode 
as part of the IRS Disks GTO program 02, AOR key 
3567360. The source was observed with three modules: 
short-hi, long-hi, and short-lo [both first and second or- 
ders] . A small raster map consisting of two positions par- 
allel to the slit and three perpendicular to the slit was 
made with each module. The raster spacings parallel to 
the slit were selected so that the usual IRS staring mode 
point source measurement - which consists of measure- 
ments at positions one third of the slit length from either 
end - could be reconstructed from the central of the three 
perpendicular positions. The spacings perpendicular to 
the slit were a little more than half of the slit width. The 
observations at each raster position were carried out with 
an integration time of 6.3 sec. 

Examination of the two dimensional spectral images 
shows the source to be compact along the slit. Thus, 
the reduction has been carried out using only the central 
perpendicular position from each module, as though the 
data were the result of an IRS staring mode observation. 
The position of this point source, as defined by the raster 
parameters, is RA = 19 h 37 m 0.89 s , Dec = +07 o 34'9.84" 



TABLE 1 
Summary of Outflow Parameters 



Dist. 

[PC] 

(1) 


Mthin 

[Mq] 
(2) 


Rco 

[PC] 

(3) 


Vco 

[km s" 1 ] 
(4) 


T CO 
(5) 


^ obs 

[M km s^ 1 yr _1 ] 
(6) 


Lobs 

[Mq cm 2 s~~ 2 yr" 1 ] 
(7) 


150 

250 


0.008 
0.02 


0.21 
0.36 


4.55 
4.55 


4.6xl0 4 
7.7xl0 4 


7.9xl0" 7 
1.3xl0~ 6 


1.8xl0 4 
3.0xl0 4 



Note. — Outflow parameters are calculated according to the procedure described in Cabrit Sz Bcrtout (1992) and arc derived from the large 
(16 ; x 12') 12 CO and 13 CO maps presented in Figures [3] and [4] 
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TABLE 2 
B335 Photometry 
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Fig. 2. — SED and spectrum of B335. Top panel: Different 
symbols represent fluxes measured at the different beam sizes 
indicated in the legend; these data are listed in Table 2. The 
magnitude of possible aperture corrections in the 3 to 10 /jm 
region can be judged by the crosses, which are photometry 
measured in a l'/2 aperture (point-source). The black line 
shows the IRS data smoothed to ~ 0.4 fim. Bottom panel: 
The grey line shows the IRS spectrum at full resolution; the 
black line shows the same data smoothed to ~0.4 /im. Vari- 
ous features are indicated in the spectrum. Both panels: The 
IRS spectra do not join because the high resolution spectrum, 
from 10 to 20 /im, has a larger slit than the low resolution 
spectrum, from 5 to 14 /im; hence the longer wavelengths 
include more extended emission. 

(all positions in this work are given in the J2000 sys- 
tem). A small amount of extended emission may have 
been missed in this process. The synthetic point source 
obser vat ions were reduced using the Cubism software 
l|Smith et alj|2007j ). The spectrum is shown in Figure [2j 
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2.2. 12 CO and 13 CO data 

The B335 region was mapped in the J=2-l transi- 
tions of 12 CO and 13 CO with the 10-m diameter Hcin- 
rich Hertz Telescope (HHT) on Mt. Graham, Arizona 
on April 22nd, 2007. The observations were made with 
the 1.3mm ALMA sideband separating receiver with a 
4 - 6 GHz IF band. The 12 CO J=2-l line at 230.538 
GHz was placed in the upper sideband and the 13 CO 
J=2-l line at 220.399 GHz in the lower sideband, with 
a small offset in frequency to ensure that the two lines 
were adequately separated in the IF band. The spec- 
trometers, one for each of the two polarizations, were 
filter banks with 256 channels of 250 KHz width and 
separation. At the observing frequencies, the spectral 
resolution was ^0.3 km s _1 and the angular resolution 
of the telescope was 32" (FWHM). The main beam effi- 
ciency, measured from planets, was 0.85 ± 0.4. 

A 16' x 12' field centered at RA = 19 ,l 37 m 00.4 ;5 , Dec 
= +07°35'29" was mapped with on-the-fly (OTF) scan- 
ning in RA at 10" s _1 , with row spacing of 10" in decli- 
nation, over a total of 72 rows. System temperatures 
were calibrat ed by the standard ambient temperature 
load method (jKutner k. Ulich|[l98ll ) after every other row 
of the map grid. Atmospheric conditions were clear and 
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Fig. 3. — Large CO J = 2 - 1 OTF channel map of B335, 16' X 12' in size, showing channels with significant emission, from 12.52 km s — 1 
to 5.04 km s — , at a ~0.3 km s — 1 resolution. Channel velocities are indicated in the top left corner and the beam size (32") is indicated 
as the filled circle in the bottom left of each panel. Contour levels are {1, 2, 3, 4, 5, 6, 8, 10, 12, 14, 16} X 0.5 K-Tt . 



stable, and the system temperatures were nearly con- 
stant, and averaged T sys = 191 K (SSB). 

Data for each CO isotopomer were processed with the 
CLASS reduction package (from the University of Greno- 
ble Astrophysics Group), by removing a linear baseline 
and convolving the data to a square grid with 10" grid 
spacing (equal to one-third the telescope beamwidth). 
The intensity scales for the two polarizations were de- 
termined from observations of standard sources made 
just before the OTF maps. The gridded spectral data 
cube s were processed with the Miriad software pack- 
age (jSault et al.lll995T ) for further analysis. This process 



yielded images with rms noise per pixel and per velocity 
channel of 0.21 K-T^ for both the 12 CO and 13 CO tran- 
sitions. The 12 CO and 13 CO channel maps are shown in 
Figures [3] and 2J respectively. 

Using the same set-up as that outlined above, we made 
a deeper map of a 2.'5 x 4.'5 region centered at RA = 
19 /t 37 m 00.8 s , Dec = +07°34'07", near the protostar. 
These observations were carried out on March 19, 2008, 
and were repeated 7 times. The average T sys = 205 K 
(SSB), and was nearly constant throughout the obser- 
vations. After combining the 14 images with weights 
proportional to the individual map rms values, we ob- 
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Fig. 4.— Large 13 CO J = 2 - 1 OTF channel map of B335, same 
size and spectral resolution as Figure [3] Channels with significant 
emission are shown, from 7.62 km s to 9.33 km s _1 . Channel 
velocities are indicated in the top left corner and the beam size 
(32") is indicated as the filled circle in the bottom left of each 
panel. Contour levels are {1,2,3,4,5,6} X 0.5 K-T^. 



tain final rms values per pixel and per velocity channel 
of cr T1 = 0.091 K ( 12 CO) and a T * = 0.084 K ( 13 CO). 
These maps are shown in Figures [5] and [H while Figure [7] 
shows position-velocity diagrams derived from the same 
data. 

3. THE DISTANCE 

A quantitative interpretation of our observations re- 
quires an estimate of the distance to B335. Here 
we examine the often-quoted and widely accepted 
iTomita et al.l (|1979l ) distance to B335 of 250 pc and up- 
d ate it with modern d ata. 

ITomita et al.l (|1979| ) derive a distance to B335 rang- 
ing from 130 t o 250 pc, ba s ed on photoelectric stellar 
data from the iBlanco et al.l (|1970f ) catalog, and a dis- 
tance of 600 pc based on star counts. They adopt 250 pc 
as a plausible value but do not provide further details or 
show figures of Ay v s. distance or of stell ar density in the 
field of interest. The IBlanco et al . (1970) catalog d ata in 
the region of B335 originates from Bouigue et all (|1961l ) 
catalog photometry. For completeness we note that the 
origin of the s pectr al types of the stars listed in the 
iBouigue et al.l (J1961I ) catalog are from an evidently un- 
published provisional spectral classification determined 
by l'Observatoire de Marseille; the quoted spectral types 
appear generally reliable but d o not always agree with 
those quoted in SIMB ALP . The IBlanco et all (|1970f) cat- 
alog references 9 stars ()Bouigue et al.lll961l) within 45' of 
B335, well outside the optical extinct i on ext ent of B335, 
r ~ 3'. Therefore, the ITomita et al.l (|1979f ) distance of 
250 pc is highly uncertain. 



We show an updated version of the Ab vs. distance 
plot in Figure [H] The data are from the on-line VizieR 9 
service. We find 32 stars within 20' of B335 with spectral 
classifications. Of these, we reject all stars without both 
B-band and K-band photometry, as well as those with 
spectral classifications of G2 or later. Later-type stars 
mostly add scatter to our analysis because a small error 
in the spectral type yields a large error in both Ab and 
distance. In all cases for which we do not have a luminos- 
ity class we assume that the star is on the main sequence. 
We a ssume A B = 1.09 x E(B - K) (|Rieke fc Lebofskvl 
Il985f ) . The res ulting data are prese nted in fig. [8] We 
also indicate the ITomita et alJ (|1979f ) distance estimates 
(at 130, 250, and 600 pc) as dotted lines. We do not 
see any noticeable feature at or near 2 50 pc. To put a 
lower bound on the distance we use the lLallement et al.l 
(2003) maps of the Local Bubble; in the direction of B335 
(7 = 44.9°, b — —6.6°) they measure the extent of the lo- 
cal low density bubble to be ~ 60 pc. We set the upper 
limit to the distance to B335 to be ~ 200 pc since in Fig- 
ure [5] we see approximately constant stellar Ab values 
outside of this distance. We adopt a distance of 150 pc 
in the following discussion, indicated as a solid line in 
Figure [5] Our adopted distance is uncertain by nearly 
a factor of two, but it falls near the middle of the plau- 
sible ran ge. It also reduces a number of discrepancies 
found bv lHarvev et al.l (|2001l ) in modeling the extinction 
of the globule, and helps reconcile the inferred properties 
for the protostar with the luminosity of the globule (see 
§5). 

4. RESOLVED STRUCTURE 

4.1. The Outflow 

The 12 CO and 13 CO maps are presented in Figures [3] 
through[7l The large maps (16' x 12') show the full spa- 
tial extent (~ 8') of the outflow as well as the velocity 
extent. The western side of the outflow is redshifted 
with respect to the line of sight while the eastern side is 
blucshiftcd. These maps are also used by Shirley et al. 
(2008, in preparation) to calculate the outflow opening 
angle. They use this angle (~55°) to exclude from their 
analysis of the submillimeter dust opacities regions where 
dust may be affected by processes such as shock heating. 
The 12 CO and 13 CO lines are well resolved with 250 KHz 
(~ 0.3 km s _1 ) channel width; in 12 CO the blucshiftcd 
lobe is dominated by velocities ranging from ~ 5.04 to 
^6.99 km s _1 while the redshifted velocities range from 
~ 9.27 to ~ 12.52 km s^ 1 (see Figure [3]). The subfield 
map (Fig. [5|) shows that the velocity extent of the red- 
and blucshiftcd emission is in fact larger; however, these 
maps do not cover the full 8' region of the outflow. The 
13 CO data do not show an outflow; the range of these ve- 
locities is much narrower than the 12 CO velocities, from 
about 9.33 to 7.62 km s _1 . The 13 CO data are domi- 
nated by colder emission: in the large map (Fig. [4]) we 
see little evidence of the outflow lobes observed in 12 CO. 
However, some east- west elongation can be observed in 
the sub-field map, see Figure O due to the lower noise 
(rms ~ 0.084 K-TJJ in the smaller map as compared to 
the full field map (rms ^0.21 K-T^). Both large maps 
show an extension of emission to the north and east of 
the protostar, on a scale of ~ 9' in 12 CO and ~ 6' in 
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Fig. 5.— 12 CO J = 2 - 1 OTF sub-field channel map of B335, 2f5 X 4.' 5 in size, centered on RA = 19 h 37 m 00.8' ! , Dec = +07°34'07'.'4, near 
the location of the protostar. Channels with significant emission are shown, from 15.44 km s -1 to 2.11 km s — 1 , observed at a resolution 
~0,3 km s _1 . The temperature scale is indicated in the right side; contour levels are {1,2,4,8,16} X 0.25 K-T^. Note the structures 
located north and south of the main east-west outflow region, visible here in the south at ~9.27 km s — 1 and in the north at ~7.31 km s — . 



13 CO, which is also observed in the 160 /iiri mosaic. The 
outflow parameters calculated from the large area maps 
are presented in Table [TJ 

The IRAC images reveal a resolved emission structure 
consisting of two outflow lobes oriented in the east-west 
direction, seen clearly in Figure Q] in all four IRAC chan- 
nels. The observation that the eastern lobe is brighter 
than its western counterpart is explained by the fact 
that the outflow is oriented slightly out of the plane of 



sky; this orientation is corroborated by our CO maps 
in which the eastern lobe is slightly blueshifted rela- 
tive to the western lobe, as discussed above. It is also 
in agreement with othe r measurements: for example, 
ICabrit fc Bertoutl (|l992t ) measured the inclination angle 
to the line of sight to be 82°. The outflow not being 
aligned exactly on the plane of the sky will give rise to 
differential extinction between the eastern and western 
lobes, causing the western lobe to be more obscured by 
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Fig. 6. — 13 CO J = 2 - 1 map; same as Figure[5] Channels with 
significant emission are shown, from 7.10 km s" 1 to 9.48 km s" 1 . 
Contour levels are {1, 2, 3, 4, 5, 6, 7, 8, 9} X 0.25 K-T^. 
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Fig. 7. — 12 CO (top) and 13 CO (bottom) position-velocity 
diagrams for the sub— field map. The 12 CO contour levels are 
{1,2,4,8,16} X 3<r, where 3<r =0.28 K-T^. The 13 CO contour 
levels are also {1, 2, 4, 8} X 3<T, where 3cr = 0.25 K-T^. The 12 CO 
north-south cut is centered on RA = 19 h 37™0.75 s ; the east-west 
cut is centered on Dec = 07°34'7'.'05. The 13 CO north-south cut 
is centered on RA = 19 h 37 m 0.83 s ; the east-west cut is centered 
on Dec = 07°34'5'/26. The vertical lines (dashed and solid arrows) 
in the 12 CO and 13 CO NS cuts indicate the offsets (80") at which 
the flattened molecular core is detected. 



the surrounding cold globule material out of which the 
protostar was formed. 

To improve the spatial resolution of the Spitzer data, 
we present the deconvolved HiRes images of B335 in 
Figure [9] (right-hand column) . We reprocessed the raw 
(BCD) data from the Spitzer archive, using a version 
of HiRes deconvolu tion developed for Spitzer images by 
iBackus et" al. ( 2005 ) based on the Ri chardson-Lucy algo- 
rithm (Richardson 1972; Lucv 1974) and the Maximum 
Correlation Method (|Aumann et al.| [l990) . HiRes decon- 
volution improves the visualization of spatial morphol- 
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Fig. 8. — A B vs. distance for stars within 20' of B335. The 
error bars indicate the magnitude of the effect of a 15% change 
in each direction in the (B — K) color. The dot ted lines indicate 
the range in distances derived by Tomita et al. (1979): 130, 250 
(adopted by those authors), and 600 pc, respectively. Note that we 
do not reproduce a discontinuity in Ab at 250 pc. The dashed lines 
indicate plausible ranges in the distance based on the measured 
extent of the local bubble and the Ag data shown here (see § 3): 
60 and 200 pc. Our adopted distance of 150 pc is indicated by the 
solid line, and represents a conservative estimate. 

ogy by enhancing resolution (to sub-arcsec levels in the 
IRAC bands) and re moving the conta minating sidelobes 
from bright sources (jVelusamv et al.ll2008l ). It achieves 
a factor of 3 enhancement over the diffraction limited 
spatial resolution from ~ 1.22A/D to ~ 0.4 A/D (e.g. 
from 2.4"to ~ 0.8"at 8/mi). Typically the diffraction 
residue in the HiRes images is at levels well below 0.1% 
to 1% for the IRAC and MIPS 24 /im bands. The Point- 
source Response Functions (PRFs), required as inputs 
to HiRes, were obtained from the SSC; we used the ex- 
tended (about 128 detector pixel size in each direction) 
versions, supplied by Tom Jarrett (private communica- 
tion, 2007). The observed images, in the form of raw 
BCDs, were first post processed to fix contaminated data 
for saturation, column pulldown, muxbleed, muxstrip- 
ing, jailbars, and outliers. In the Richardson-Lucy algo- 
rithm the best resolution enhancement is attained when 
the background is zero. Accordingly, we used the back- 
ground subtracted BCDs. In HiRes we used 50 iterations 
for the IRAC bands and 100 for MIPS, which resulted 
in a FWHM of 0.55"to 0.75"for IRAC channels 1-4 and 
1.6"at 24 /im for a point source in the deconvolved im- 
ages. All bands are relatively clean from artifacts except 
IRAC channel 4 (at 8/im) which has some uncorrected 
muxbleed. In the HiRes images diffraction residues near 
the Airy rings are low (at a level < 0.05% of the respec- 
tive peak intensities at 3.6, 4.5, and 5.8 /im) providing 
enhanced visualization closer to the protostar. 

In the HiRes 5.8 /im image we measure the outflow 
opening angle to be about 70°, while in the deconvolved 
8.0 /im image we measure the angle to be between 62° 
and 67°. Shirley et al. 2008 (in preparation) measure 
the 12 CO outflow opening angle to be closer to 55°. In 
Figure [TU] we show the 12 CO integrated contours over- 
laid on the IRAC 4.5 /im minus 3.6 /im image. It ap- 
pears that the molecular outflow is more confined than 
if it were simply expanding from the opening angle in- 
ferred from the IRAC images. Also note the presence of 
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Fig. 9.— Gallery of images of B335 (CB199) at the indicated 
wavelengths. The images in the right column have been processed 
with HiRes for improved spatial resolution. North is up and east 
is to the left. 

Herbig Haro (HH) objects in the subtracted image along 
the inner walls of the confining cavity, where it is proba- 
ble that dense shocked gas produces molecular hydrogen 
and CO emission that accounts for much of the signal in 
the IRAC images. These emission features fall in IRAC 
channel 2. The HH objects (some having a striking fil- 
amentary geometry) imply that the inner cavity wall is 
quite sharp. More detailed discussion of the B335 HH 
objects can be found in, e.g. JGalfalk fc Olofssonl (|2007j ) 
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Fig. 10. — IRAC 4.5 /im minus 3.5 /im image with 12 CO inte- 
grated contours at the indicated temperatures. The countours are 
integrated over to 15 km s — 1 . The beam size (32") is indicated 
by the grey circle. Note the presence of HH objects, knots of dense 
shocked material producing molecular-H emission, located near the 
southern and eastern regions of the 4 K km/s contour. North is up 
and east is to the left. 



and lGalvan-Madrid et all (J2004D . 

The IRS spectrum (see Figure [2]) shows emission lines 
at ~ 26 /im and ~ 35 /im which we identify as [Fell] 
and [Sill], respectively. These lines may be produced in 
the interaction between the outflow and the surrounding 
neutral material. We note that they are also seen in the 
spectrum of the north-e ast lobe (HH47A) in the HH4 6/47 
outflow, as discussed b y lNoriega-Crespo et all (|2004f ) and 
IVelusamv et al.l (|2007l ) . We also observe other features of 
interest in the IRS spectrum, indicated in Figure^ like 
H 2 Q and CO2 ices a nd a deep silicate absorption feature 
(jvan Dishoeckll2004ft . 

4.2. The Shadow 

The IRAC images show a shadow, located near the out- 
flow lobes. In Figure Q] the shadow is evident at 3.6 /im 
and 8.0 /im, just to the south of the outflow lobes, as a 
dark nearly circular feature about 20" in extent. We also 
marginally detect a northern counterpart in the 3.6 /im 
and 8.0 /im images (see Figures fTI and ITT]) . We show the 
IRAC images on a log scale, with the minimum and max- 
imum image values set to highlight the shadow and the 
outflow lobes; no other special processing of the images 
was used to display the images. 

The likely cause of the IRAC shadow is dense globule 
(dust and molecular) material near the protostar. The 
observed difference in brightness between the east and 
west outflow lobes (apparent at all IRAC wavelengths, 
see Figures [TJ and |H|) , is explained by: (a) the outflow 
orientation with the eastern lobe pointing slightly to- 
wards the observer, and (b) the presence of dense ma- 
terial at large radii (compared to the ~ 200 AU circum- 
stellar disk) causing differential extinction between the 
two lobes. This interpretation is supported by the spa- 
tial ^disWbution and_measured color excesses observed 
by lHarvev et al.1 ((2001). These data are presented in 
Figure Q2J where we show the locations of the stellar 
sources (circles) overlaid on the binned 8.0 /im image. 
Here t he colors indicate a rough extinction distribution 
for the I Harvev et al.l (|2001h sources: red circles indicate 
sources with (H — K) color excess greater that 2.9 (equal 
to the mean value plus 2x the rms in the color excess 
distribution), the orange circles indicate color excesses 
between 2.1 (equal to the mean value plus Ix the rms 
in the color excess distribution) and 2.9, and the yellow 
circles indicate sources with excesses less than 2.1. In 
this figure, the ellipse (with a 3:2 major to minor axis 
ratio) indicates the rough shape of the hole devoid of 
stars near the center of the cloud. The location of the 
most obscured sources around the outer circumference of 
the southern shadow is suggestive that the shadow is in 
fact tracing a flattened overdensity in molecular material. 
On the northern side the shadow is also detected in the 
IRAC image; its location, like the southern counterpart, 
is bounded by background stellar detections, indicating 
that the molecular core is in fact flattened and that it 
has a significant degree of symmetry. 

In the text that follows we analyze the shadow and 
derive optical depth and density pr ofiles; we use tech - 
niques similar to those presented in lStutz et al.1 (|2007l) . 
The shadow is roughly centered on RA = 19 h 37 m 01 s , 
Dec = +07°33'55". Although the shadow is observed at 
both 3.6 /im and 8.0 /im we apply the following analysis 
only to the 8.0 /im shadow because the 3.6 /im image 
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Fig. 11. — IRAC 8.0 fira image of B335, with inset, both shown on a log scale and at the original mosaic pixel scale of l'/2 pix . The 
8.0 /im shadow can be observed near the center and just south of the protostar. We also show the HiRes 8.0 (im deconvolved contours 
in red (west) and blue (east); these colors correspond to the CO outflow kinematics. The orange ellipse has a 3:1 major to minor aspect 
ratio and indicates the rough size and shape of the flatte ned molecular cor e. Inset: The solid black contours indicate the location of the 
circumstellar disk (~ 200 AU) observed at 1.2 mm by Harvey et al. (2003a); the black dashed circle (with 3" radius) indicates the 24 fj,m 
PSF. North is up and east is to the left. 



has too many complicating factors, namely contamina- 
tion from the nearby bright source just west of B335 and 
many more stellar sources in and around the region of 
interest. In summary, to derive an 8.0 /im profile we go 
through the following steps: (1) re-normalize the image 
to our best estimate of the true background value (Fig- 
ure [13]), (2) divide the image up into regions inside and 
outside of the shadow (Figure [H)) . (3) use these regions 
to derive optical depth, density, and extinction profiles, 
all of which are equivalent, the latter two being subject 
to an assumption about the dust properties at the wave- 



length of interest (Figure [TS)) . and (4), compare the de- 
rived pr ofile to available ste llar color excesses, in this 
case the lHarvev et al.l (|2001l ) data, see Figure [HJ Here 
we discuss this process in more detail. 

The background level in the images, which we term 
foe an d which is composed of instrumental background 
and zodiacal background, must be measured and sub- 
tracted from the images to analyze the shadow. We esti- 
mate foe by adjusting the level so that the depth of the 
shadow and the extinction inferred from stellar observa- 
tions agree. We compare this level with the pixel values 
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Fig. 12.— IRAC 8.0 /im image of B335, binned 2x2 for a 
2''4 pix - 1 resolution. The 8.0 /im shadow can be observed near 
the center and just south of the outflow lobes; a northern shadow 
counter part is marginally de ctected. The circles indicate the loca- 
tions of Harvey ct al. (2001) near-IR stellar sources: the large red 
circles indicate sources with (H — K) color excesses greater than 
2.9 (2cr, where for this set of observations cr = 0.8), the medium 
orange circles indicate sources with color excesses between 2.1 (lc) 
and 2.9, and the small yellow circles indicate sources with excesses 
below 2.1. The cross indicates the location of the Harvey ct al. 
(2003a) detection of the central source. The dashed ellipse (shown 
with a 3:2 axis ratio) is meant to highlight the shape of the inner 
region without stellar detections. North is up and east is to the 
left. 
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Fig. 13. — IRAC 8.0 /im pixel value distribution for a region 
2'2 on a side, 2'2 north of B335. The solid line indicates the 6th 
percentile value that we use for the image background value, equal 
to 827.408 MJy sr~ x , chos en to get the b est match to the stellar 
extinction values from Harvcv ct al. (2001). 

in a 2f2 squared box located 2f2 north of the source se- 
lected for minimum surface brightness; the distribution 
of pixel values is shown in Figure fT3l The chosen value for 
Jdc = 827.40 MJy sr _1 is the 6th percentile pixel value 
in the box. The reasonably good agreement between our 
detected background and the assumed one validates our 
procedure. 
We divide the image into nested and concentric ellipti- 




FlG. 14. — IRAC 8.0 /an image of B335, with two stars masked 
out (black boxes). The blue ellipses, centered on the Harvey et al. 
location for the protostar (marked with a green circle) indicate the 
nested and adjacent regions we define to analyze the shadow. The 
green lines indicate a 60° wedge between position angles 240° and 
300°, outside of which we reject all pixels due to light contamina- 
tion from (a) the protostar reflection lobes, and (b) the diffraction 
spikes from the nearby very bright source. 

cal regions sho wn in Figure 1 1 4 1 The elliptical regions are 
centered on the lHarvev et al.l (|2003a[ ) location of the pro- 
tostar, and have a major-to-minor axis ratio of 3:1. The 
ratio was chosen to match roughly the location and shape 
of the 8.0 /im shadow and to satisfy the mininum require- 
ments of the observed geometry: inspection of Figure [141 
provides the constraint that the flattened core must lie 
in front of the entire western outflow lobe while still pro- 
viding the obscuration observed as a shadow at 8.0 /im. 
Smaller axis ratios (rounder ellipses) will violate these 
constraints. Similarly, larger axis ratios may allow for 
some of the emission from the outer region of the western 
lobe to appear unobscured. We reject all pixels outside 
of a 60° wide wedge pointed directly south; this wedge is 
indicated as two green lines in FigurcQj] It was chosen to 
avoid the diffraction spikes of the nearby bright star and 
to mask out regions affected by emission from the outflow 
lobes. As can be seen in Figure[TH we also mask out two 
bright sources (black boxes) that would have affected the 
shadow profile. We then tabulate the mean pixel value in 
each elliptical region, plotted in Figure [T5l (upper panel, 
square symbols). We use this distribution to calculate 
an estimate of the optical depth, r = — ln(//io), hi the 
shadow. To do so we must measure Io, t ne unobscured 
flux level near the shadow. We measure Iq from the re- 
gions where the mean pixel values approach a constant, 
in this case at elliptical major axis values greater than 
56". We average over all pixel values outside this region 
and inside 78"; we measure the value of Iq = 8.26, in- 
dicated in the upper panel of Figure [15] as a solid line. 
In the lower panel of Figure [TBI we plot exp(— r 8 ) = I / Iq 
versus major axis distance. One can clearly see that the 
inner-most region is greatly affected by emission from the 
outflow lobes. 

Motivated by the concern that emission from the bright 
source just to the west of B335 could be affecting our 
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Therefore the nearby bright source is not significantly 
affecting our shadow profile. 

Finally, using this r-profile, we derive a n extinction 
profile. Following the convention outlined in lStutz et al.l 
(|2007f ). and assuming an 8.0 (im dust opacity of n a bs,s = 
8.13 x 10 2 cm 2 gm" 1 calculated bv lDraind (|2003al |bh from 
his Ry = 5.5 model, we derive the Ay profile plotted in 
Figure 1161 Here, the average mass column density in 
each elliptical annulus is given by 



~R 



Kabs,8 



■/• 



(1) 



where /(= 100) is the gas-to-dust ratio, and K a bs,a is the 
absorption cross-section per mass of dust, and has the 
assumed value noted above. The choice of a model with 
a high Ry value relative to the diffuse ISM is intended 
to account for the grain coagulation expected in dense, 
cold cores. The extinction profile is given by 
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Fig. 15. — IRAC 8.0pm shadow profile. Top panel: Aver- 
age pixel values as a function of major axis distance from the 
lHarvev et a l. (2003a) location of the protostar. Solid, short- 
dashed, and dotted lines indicate, respectively, the derived 
unobscured flux levels (7o) in the three indicated regions: the 
full 60° wedge, the eastern 30° half, and the western 30° half 
(located at the position angles indicated in teh top panel). 
The vertical long-dashed line indicates the six outer regions 
used to derive Iq. Bottom panel: exp( — r) = I/Ip as a func - 
tion of major axis distance from the lHarvev et al.l lj2003al ) 
location of the protostar. Both panels: Squares indicate the 
profile derived using the full 60° wedge, triangles indicate the 
eastern 30° half of the wedge, and crosses indicate the west- 
ern 30° half of the wedge. The bounding position angles of 
these wedges (240°, 270°, and 300°) are indicated in the top 
panel. Note that there are not significant differences between 
the profiles, indicating that the derived profile is not strongly 
affected by the nearby bright source. 



derived shadow profile we perform the following test. We 
divide the 60° wedge between position angles 240° and 
300° into two 30° wedges, one on the east side, away from 
the bright source, and one on the west side. We proceed 
with the same analysis as described above and plot the 
resulting profiles in Figure [T5l These profiles (triangles 
indicate the east wedge, crosses indicate the west wedge) 
do not show significant variations between each other, or 
relative to the full wedge (again, indicated as squares). 



We have assumed the conversion factor from the re- 
lation N(H 2 )/A V = 1.87 x 10 21 atoms cm" 2 mag -1 
(Bohlin et al. 1978), where we assume that all of the hy- 
drogen is in molecular form. We have adopted Ry = 3.1 
— a value appropriate for the diffuse ISM where th i s rela - 
tion was measured — to convert the iBohlin et al.l (|1978T ) 
relation from E[B — V) to Ay. At longer wavelengths 
than V the extinction law does not appear to be a strong 
function of density. Here /j,h 2 = 2.8 is the effective molec- 
ular weight per hydrogen molecule. For this set of con- 
stants Ay — 14.04 x Tg. Following a different method, 
where t he conversion is obtained directly from the lDrainei 
(|2003af bT) dust models, Ay = K a bs,v / K a bs,s x ^8 = 
9.7 x rg. The two Ay determinations are consistent 
within the uncertainties of the measured infrared ex- 
tinction values (see below). To take possible freezeout 
and grain growth effects int o account we consider the 
lOssenkopf fc Henningj (|1994[ ) model opacities; for mod- 
els generated at a density of n = 10 6 cm -3 these authors 
derive 8.0 /mi opacities of K a bs,a = 1.06 x 10 3 cm 2 gm" 1 
and K a bs,$ — 1-21 x 10 3 cm 2 gm" 1 for thin and thick 
grains, respectively. These opacities will diminish the 
derived Ay values by factors of 1.3 and 1.5 respectively. 
We n ote that the range in m easured infrared extinction 
(e.g.. lIndebetouw et al. 2005; Fl aherty et a l. 2007) is sig- 
nificant; the Ay determination is therefore uncertain by 
a factor ~ 2 and as a result so is the column density. 

Figure [T|)] shows the derived shadow extinction profile 
(square symbols ); we co mpare our profile to the values 
derived from the lHarvev et a l. (2001) stellar data in each 
elliptical region. To make use of the st ellar observed 
H — K) colors we follow the procedure in lHarvev et al.l 
2001). In addition to the statistical uncertainty in 
(H — AT), they also measured the uncertainty due to the 
intrinsic scatter in (iJ — K) of background stars deter- 
mined from a nearby "OFF" field. The mean intrinsic 
(H - K) = 0.13 ± 0.16 for the "OFF" field. To deter- 
mine the color excess, the stellar (H — K) colors have 
0.13 subtracted from them; to determine the total errors 
we add 0.16 in quadrature to the statistical errors. Fi- 
nally, we use Ay/E(H — K) = 15.87 to convert to Ay 
(|Rieke fc Lebof skv 1985). To make a fair comparison be- 
tween the shadow profile and the stellar color exess, we 
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Fig. 16. — Ay values fo r the 8.0 /im shadow (squares) compared 
to the lHarvev et a l. (2001) stellar Ay values (circles) for each ellip- 
tical region. We use the six outer radii — indicated by the dashed 
line — to normalize both the shadow and stellar extinction profiles. 



must subtract off the mean value of the stellar extinc- 
tion in the regions used to find Iq for the shadow, i.e. we 
account for the differences in the two zero-points for the 
profiles. In this case, (A v ) = 14.52, or E(H-K) = 0.92, 
for the stars in this region. 

The maximum extinction level derived for the shadow 
is Ay ~ 100. Our spectrum (Fig. [2]) shows silicate ab- 
sorption with a r > 2.5. Because of possible complexi- 
ties with interpreting this spectr um, we consider r ~ 2.5 . 
Applying a standard conversion (jRieke fc Lebofskvl 19851 ) 
the corresponding extinction is Ay > 40. We interpret 
the measurement as a lower limit because radiative trans- 
fer in a disk with a temperature gradient will tend to fill 
in the intrinsic feature. If we assume an equal amount of 
material behind the star, the limit on the total column is 
equivalent to Ay > 80. Thus, the spectrum is consistent 
with our shadow analysis: there is a column density of 
N(H2)~ 10 23 cm~ 2 or more towards the protostar. For 
a radius of 20", this column yields a mass ~0.5 M©. 

Figure [11] shows additiona l components of t he source 
in relation to the shadow. lHarvev et al.l (|2003af ) used the 
IRAM interferometer in the A configuration, which has 
baselines of 43 to 400 m. At this high resolution they 
detect the protoplanetary disk in emission, with a max- 
imum diameter of about 1" (150 AU, see Fig. [TTj) . but 
resolve out larger structures. lHarvev et a l. (2003b) used 
the D configuration with baselines of 15 to 80 m, which 
allow for larger structures to be resolved. Their best fit 
model to the data (p — 1.65, a = 40°, see their Figure 3) 
shows a structure about 7" in extent aligned in the north- 
south direction (not shown in Fig. Ill[) . This emission 
likely originates from the innermost material associated 
with the 8.0 /im shadow. The dashed circle indicates the 
24 /mi PSF (~ 3" in radius) which is coincident with the 
protostar and disk. 

Our analysis does not include any treatment of the pos- 
sible foreground emission that may arise from the surface 
of the globule its self; this emission, if present, would 
have the net effect of causing an underestimate of the 
derived optical depth. A rigorous modeling of this emis- 
sion component is beyond the scope of this paper. How- 
ever, inspection of Figure [16] indicates that if this effect 



is present it is not large; the lHarvev et al.l (|2001l ) stellar 
data lie preferentially above our shadow profile, at higher 
Ay's. Foreground emission could explain this systematic 
offset. However, other uncertainties, discussed above, are 
likely to dominate our extinction profile estimate. 

4.3. Detection of a flattened molecular core 

The sub-field 12 CO and 13 CO channel maps show a 
north-south structure that appears to be rotating, seen 
in Figures [5] and [6] This structure extends out to a ra- 
dius of about 80" (12000 AU) from the protostar. In 
Figure [5] the southern component is detected between 
roughly 8.94 km s^ 1 and 9.92 km s _1 ; at 8.94 km s _1 the 
image is dominated by the center velocity emission. The 
northern counterpart is weaker in emission, but is appar- 
ent at 7.31 km s _1 and fades at velocities smaller than 
6.66 km s . The maximum relative velocity between the 
two sides is therefore about 3 km s -1 . The 13 CO map 
shows little indication of this velocity gradient; however, 
there is a north-south elongation of the emitting region. 
Additionally, CO freeze-out may complicate the analysis 
of the molecular data. The difference in appearance in 
12 CO between the northern and southern regions can be 
explained by clumpy material in a rotating structure, ex- 
pected by models (see below). Thus, both our extinction 
shadow image and our CO line mapping indicate that the 
dense core of the globule is flattened and may be rotating 
around the protostar. 

There have been many indications of this flattened 
structure previously, on different angular scales depend- 
ing on the type of observat i on. Fo r example, the 800 /mi 
image by iChandler et al.l |l990) indicates a low sur- 
face brightness north-south extension of about 2', com- 
pared with only abou t 1' for the east- west extension. 
IVelusamv et al.l ()1995l ) find a very similar structure in 
their filled aperture image in CCS at 22 GHz. They also 
present channel ma ps that show no convincing evidence 
of overall rotation. iHodappI (|1998f ) shows deep images 
at H and K that reveal shadows of the cloud against the 
diffuse light from the background, a technique analogous 
to ours at 8.0 /im. At H, the shadow is slightly flat- 
tened, but it is significantly more so at K-band. This 
change is probably due to the ability of the K-band im- 
age to penetrate larger optical depths than at H-band. 
ISaito et all (]1999h show a similar structure in C 18 and 
find marginal evidence that it is rotating with a small 
redshift to the north and blue shift to the south. The in- 
dicated rotation is only about 0.05 km/sec over the range 
o f 2' (0.09pc). 

IWilner et all (|2000D observed in CS J = 5-4 with the 
IRAM interferometer and an angular resolution of 2.5"; 
their channel maps do not appear to show the full struc- 
ture, but there is a signal near the systemic velocity that 
is well resolved with a length of 5" to 6" oriented ap- 
proximately N-S. This structure is also seen in the IRAM 
interferom eter continuum data o btained at 1.2mm simul- 
taneou sly (jHarvev et al.ll20d3bh and by |j0rgensen et all 
(2007) with the SMA (the latter image also shows a 
faint extension to the north and another to the SSW). 
It is plausible that the latter three references underesti- 
mate the extent of the larger, low surface bright struc- 
ture because it is partially resolved out with the inter- 
ferometers. This possibility is partially supported by 
the strong detection of the north-south structure over 
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TABLE 3 

Mass Estimates 



Data 



(1) 



Temp. 

[K] 

(2) 



Radius 



(3) 



Mass [M Q ] 

[M ] 

(4) 



8.0 jum 

160 (im 
160 (im 
160 (im 
160 (im 
160 fira 
160 (im 
13 CO 

13 co 

12 C0 6 



10 
10 
12 
12 

11 
11 
30 
10 



20 
18 
80 
18 
80 
18 
80 
300 
300 
80 



0.5 
13 
69 
10 

15 
3.2 
5.2 
^0.08 a 
^0.2 a 
35 



a Possible CO frcczcout would imply that this value is a lower limit. 
Enclosed mass, calculated assuming a circular velocity — 1.63 km s — 
(see §4.1) 

a l ength scale of 2' in the filled aperture CS observation 
bv lVelusamv et al.l (|1995l ). 

Three-dimensional numerical simulations of cloud col- 
lapse are growing rapidly in sophistication. In general, 
they show that the collapse proceeds quickly into fila- 
ments, or disks with spiral arms, and with typical di- 
mensions of 1000 AU (diameter) flW hitcho use fc Bate! 
l2006HKrumholz et al.ll2007t I Arreaga-Garcia et al.ll2007l) , 
a res ult that was also antici pated in some earlier work 
(e.g.. lNelson fe Langerjll997l ). In the simulations, these 
structures often break up into binary or multiple form- 
ing stellar systems with dense and more compact cir- 
cumstellar/protoplanetary disks. In any case, the in- 
terstellar medium becomes organized into two distinct 
classes of structure: 1.) an evolving, pseudo-stable pro- 
toplanetary disk of size up to a couple of hundred AU; 
and 2.) what we term a flattened molecular core to in- 
clude the variety of filaments, disks, and spiral arms 
produced by the simulations, and of typical size 1000 
AU. Viewed edge-on, most versions of flattened molec- 
ular cores would appear as linear structures similar to 
edge-on disks. By an age of order 100,000 years or less, 
the natal cloud dissipates to reveal the new protostar, 
and the flattened molecular cores have disappeared. The 
protostell ar disk remains, within a ty pical 200 AU diam- 
eter (see lAndrews fc William s 2007). In turn, it evolves 
and dissipates on a time scale of a few million years, 
leaving a planetary system with typical dimensions of 
200 AU. 

The B335 structure would appear to be similar to the 
elongated s hadow in the IRAC images of L 1157 dis- 
covered by lLoonev et al.l ((2007). However, the L 1157 
shadow is ~ 15,000 to 30,000 AU in diameter , some what 
larger than that around B335. lLoonev et al.l (2007) sug- 
gest that it represents a flattened density enhancement 
produced as a result of the general influence of magnetic 
fields or rotation in the collapse of a cloud; that is, it is 
part of the general trend that yields a flattened molecular 
core, but probably on a size scale that decouples much 
of its volume from direct participation in the formation 
of the central protostar. 

4.4. Globule Mass 

A number of lines of evidence suggest that the full 
globule has substantial mass. We find ~ 0.5 M Q within 



a 40" beam, bas ed on the absorbing column of material. 
iGee et ail (J1985I ) modeled the spectral energy distribu- 
tion (fitting it with a 14 K temperature and a v 2 emis- 
sivity) to find a mass that translates to ~ 4 Mq at our 
preferre d distance of 150 p c. Similar masses were de- 
rived bv lKeene et all (|1983 ) and iDavidsonl (| 1987ft . while 
somewhat smaller value s (by facto rs of ~ 3) have bee n 
derived bv lWalker et~all (fl990l ) and lShirlev etaTl (|2r502ft . 
Correcting these measurements to a scale comparable to 
the 12 CO source diameter (e.g., according to the relative 
flux density versus aperture at 160/iin) increases them 
by a factor of 2 to 3. 

At the same time, the calculations depend critically on 
the assumed temperature of the dust. An upp er limit of 
~ 14K can be adopted on the basis of the SEE) (|Gee et al.l 
11985ft . plus expectati ons for dust temper atures due to 
heating by the ISRF (|de Luca et al.lll993ft . In Table 3, 
we show the dependence of calculated mass on dust tem- 
perature within this constraint, both for a 40" aperture 
and for the full source. We conclude that the mass of the 
entire globule is likely to be greater than 5 M Q , perhaps 
even by a factor of a few. 

A total cloud mass of 10 M Q , if rotationally supported 
against collapse, would have a net velocity gradient of 
~ 1.5 km s . Different molecular lines give differing 
measures of the velocity gradient in the cloud, from zero 
to ~ 3 km s _1 (the latter from our 12 CO data). From 
these discrepancies it is likely that the cloud is not simply 
in smooth rotation but probably has a compl ex structure 
in ter ms of both density and dynamics (e.g. IZhou et al.l 
1990). Nonetheless, the overall gradient we observe in 
12 CO supports the argument for a mass of ~ 10 M Q or 



more. 



5. THE PROTOSTAR 



Figure Q] shows our three MIPS B335 mosaics. Fig- 
ure[9]shows the deconvolved HiRes 24 and 70 /xm images. 
In Fi gure \\A\ we indicate the location of the protostar 
from lHarvev et al.l (|2003af ) with black contours (RA = 
19 ft 37 m 0.89 s , Dec = +07 o 34'10.9"). The dashed circle 
indicates the 24 /mi PSF size, which is well centered on 
the radio continuum measurement. The 70 /jm PSF is 
well centered on the protostar, while the 160 /xm point- 
source is slightly offset most likely due to artifacts in 
the mosa ic. The SCUBA 45 and 850 /xm observations 
of B335 (Shir ley et al.l I2000D are also well centered on 
the radio observations. This object thus coincides with 
a bright point source seen in each of the MIPS bands. 
As illustrated by Figure [TT1 the protostar is invisible in 
IRAC channel 4 (and all other IRAC channels). How- 
ever, in the MIPS bands it dominates the emission. We 
detect a small asymmetry in the 24 /xm PSF in the east- 
west direction. This asymmetry becomes more evident in 
the HiRes deconvolved images. The elongation is clearly 
seen along the east-west direction after a PSF subtrac- 
tion to take out the bright symmetric core of the PSF. 
We note that the direction of the 24 /xm PSF asymetry 
(PA = ~ 100°) is somewhat misaligned with the IRAC 
emission (see Fig. |9|) but in agreement with the 12 CO 
outflow orientation (see Fig. IT0|) . The protostar accounts 
for almost all of the 70 /xm radiation and for the point 
source seen at 160 /xm as well. As noted previously, the 
160 /xm emission is also extended to the north-west of 
the YSO. This extension is in the direction of the main 
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FlG. 17.— Red Digital Sky Survey (DSS) image with 160 fim 
contours overlaid in red. The DSS image is displayed on a linear 
scale, centered on RA = lg^^-OO^ 3 , Dec = +07°36'10.0" (2' 
north of the protostar). Contour levels are {0.2, 0.5, 0.8, 1.1, 10} X 
1000 m Jy/arcsec 2 . 

body of globule, and can be seen in our 12 CO and 13 CO 
maps (see Figures [3] and HI . 

By combining the measurements across the Spitzer 
bands with those at longer wavelengths from the Kuipcr 
Airborne Observatory (KAO) and ground-based tele- 
scopes, we have produced the composite spectral energy 
distribution for B335 shown in Figured The Spitzer im- 
ages show that the infrared emission from B335 is con- 
fined spatially out to 70 /im. At 160 ^rni, the images 
show a compact source plus a diffuse component; these 
results are consistent with those of Keene (1983) show- 
ing that the flux from the source starts increasing with 
measurement aperture only longward of 140 /im. 

In Table we list the present photometric results as 
well as selected results taken from the literature. IRAC 
data are tabulated at both 2'.'4 and 24" aperture diame- 
ters to facilitate the model comparison discussed below. 
At 160 /im and beyond we also give flux densities in two 
apertures. The flux density at 160 /im is 68 Jy in a 50" 
beam and 156 Jy into a 180" beam. The flux measured 
at 70 and 160 /im is consistent with in the uncertainties 
with the earlier results published bv lKeene et alj (|1983l ) 
in this wavelength range based on measurements from 
the KAO. The flux density measurements at 350 /im and 
beyond are measured at 40" beams; measurements with a 
180" beam are not available at these submillimeter wave- 
lengths. To obtain a self-consistent estimate of the total 
flux radiated in a 180" beam, we scale the flux measured 
at 350 /im and beyond into a 40" aperture by a factor 
of 2.3, which is the factor by which the 160 /im flux in- 
creases between the two aperture sizes. This procedure 
does not contradict any of the larger beam measurements 
beyond 350 /im. In addition, as the luminosity in a v¥ v 
sense is dominated by the 160 /im luminosity, our esti- 
mates will be valid. 

Integrating over the SED plotted in Figure [2] for the 



point source component, we find a luminosity of 3.2 L at 
an assumed distance of 250 pc, consistent with previous 
results. Adopting a distance of 150 pc, the luminosity 
becomes 1.2 Lq. The larger beam data give a luminos- 
ity over a 180" region of 1.9 L , peaking at 160 /im, for 
a distance of 150 pc. The difference, 0.7 L , is consis- 
tent with heating by the interstellar radi ation field: the 
interst ellar radiation density is given by ICox h Mezgerl 
(1989) as 0.5 eV cm~ 3 . The interstellar flux incident 
on a sphere of radius ~ 0.07 pc (corresponding to 90" 
at a distance of 150 pc) is, coincidentally, 0.8 L©. The 
opacity of B335 is such that we can expect most of this 
flux to be absorbed and reradiated in the far infrared, 
thus producing diffuse emission from the cloud at a level 
consistent with our observations (see Fig. [T7| . 

Detailed modeling of the present data set is deferred 
to a subsequent work. However, to place B335 into 
context, we have compared its properties with those of 
the grid of models r ecent ly ma de available on the we b 
by IRobitaille et~aTI (|2006h and iRobitaille et~aTl (|2007h . 
These models for young protostars include disk and en- 
velope components, scattered light, and outflow driven 
cavities. The SED is computed for a range of viewing 
angles, and the distance to the source can be varied as 
well. 

We input into the SED modeler the stellar mass range 
of 0.3 to 0.5 M Q . The range of selected models is plausi- 
ble and the general shape of the model SED's match the 
observed photometry. We also fit the observed photom- 
etry (shown in Table 2) with a range of distances from 
100 pc to 300 pc. The fitter prefers models with lower 
distances and high inclinations, which corroborates our 
distance estimate and agrees with the observed geometry 
of the source. The best-fit model has a % 2 = 561, a mass 
of 0.3 Mq, a distance of 100 pc, an inclination ~ 87°, 
and a di sk outer radius of 14 5 AU, in good agreement 
with the lHarvev et al.l (|2003aT l measurement. If B335 is 
placed at the traditional distance of 250 pc, the calcu- 
lated luminosity exceeds that expected from a ~0.5 M 
protostar and it is difficult to find a consistent solution 
without invoking asymmetric output of the far infrared 
flux. The set of 10 models with lowest x 2 values generally 
agree with this best-fit scenario, with stellar masses rang- 
ing from ^0.2 to 0.7 M . Masses derrived from fitting 
SEDs to models are currently very uncertain; however, 
these masses fall in an interesting range and are only 
~5% of the total globule mass. 

6. CONCLUSIONS 

We report and analyze new Spitzer 3.6 to 160 /im data 
observations of B335 along with new ground based data 
in 12 CO and 13 CO. 

1.) We revise the estimated distance from 250 pc to 
150 pc. Our preferred di stance is more comp atible with 
the extinction models of lHarvev et a l. (2001 ) and with 
mode ls of the embedded protostar (|Robitaille et al. 2006, 
1200% . 

2.) We discover an 8.0 /im shadow and derive a density 
profile for radii ranging from ~3,000 to 7,500 AU. The 
shadow opacity at maximum is equivalent to Ay = 100 
or more. 

3.) We detect a rotating structure in 12 CO on scales of 
~ 10000 AU which we term a flattened molecular core. 
This structure is coaxial with the central circumstellar 
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disk. 

4.) It is very likely that the 8.0/xm shadow and the flat- 
tened molecular core originate from the same structure. 
5.) We measure the luminosity of the protostar to be 
1.2 L Q (at 150 pc),in good agreement with models for 
young stars of the expected mass (<~0.4 M Q ) and general 
characteristics of the object. 

As the prototypical isolated, star-forming dark globule, 
there is a huge literature on Barnard 335. Nonetheless, 
our mid- and far-infrared data from Spitzer and CO ob- 
servations from the HHT provide new insights regarding 
the configuration of the collapse of the globule and the 
properties of the protostar. We have found that the new 
observations are roughly consistent with recent models. 
Further comparisons with the full suite of observations 
have the potential to improve our understanding of iso- 
lated star formation significantly. 
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